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(54) An active vibration control method and apparatus 

(57) An Improved active vibration control system 
(100) using feedback and pseudo-feedf onward sensor 
Inputs is provided for solving the problem of random and 
repetitive active vibration control and noise cancellation 
in a system. An artificial neural network (1 1 0) is used for 
leaming the dynamics of a structure (140) and for pro- 
viding output signals that follow the state variables of the 
structure. The neural network includes a plurality of neu- 
rons which obtain biasing inputs derived from inputs to 
sensors (144, 149), as well as inputs from the other neu- 
rons in the network. Further, each neuron obtains a feed- 
back input from itself. Each Input to a neuron is weighted 
using a weighting function derived on-line. The neural 
network (110) supplies structure parameters and state 
variables to an optimal controller (1 20) which derives and 
provides a control signal to the actuators (114) so as to 
counteract vibrations and/or noise sensed in the system. 
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Description 

The present invention reiates to active vibration and noise control systenns. More particularly, tlie present invention 
relates to tlie active vibration control, wherein the control signal sent to the actuator(8) is optimized, on-line. 

5 

Description of the Prior Art: 

Vibration and noise control, although an old problem, has become an increasingly important issue in recent times, 
especially in the industrial environment This is mainly due to the increase in usage of machinery in virtually all aspects 
10 of our lives, especially with the recent introduction of office automation. Today, one majorfaclor in evaluating any product 
is its level of vibration and noise generation. Space platforms, electronic computer equipment, high precision manufac- 
turing machinery, portable generators, appliances, automobiles, medical equipment, etc. are all examples of structures 
that may experience vibration and noise generation problems. 

Formerly, vibrations were controlled using passive vibration control systems. However, passive vibratton control 
IS mechanisms suffer in that they are Incapable of handling variable speed drive systems and random excitations, especially 
at low frequencies. There have been a number of active control techniques developed. Many have been adapted to 
. compensate for the problem of noise control, while only a few have addressed the issue of vibration control. 

Much of the former work In the active vibration control area can be divided into three rr^in groups: (1 ) those that 
require a synchronizing or "synch" signal and use adaptive filtering techniques to synthesize the wave; (2) those that 
20 require direct measurement of the excitation source and may or may not use filtering to compensate for other dynamics 
such as the actuators; and (3) those that perform off-line design analysis and assume, erroneously In many cases, that 
signals and systems do not change with time. 

One problem with systems requlririg a synch signal is that those systems have no mechanism to compensate for 
nonrepetitive (random) vibrations or noise sources. As to the second group above, systems which require direct meas- 
2S urement of the excitation (a direct feedforward signal) may not be possible or nnay be impractical to implement, due to 
the difficulty in locating a sensor at an excitation point. Further, equipment useful for direct measurement of an excitation 
source may potentially increase the cost of the system. 

Further, one problem inherent with any system that requires off-line modeling and design is that real systems change 
over time, especially If Inputs to the system are time varying. A system designed off-line cannot take changes in the 
30 system Into account. As such, none of the above-listed approaches provides an economical, practteal and efficient 
solution to the problem of vibrations in a system. 

Additionally, many prior art systems are slow due to the time Intensive calculations used to obtain a transfer function 
used by the controller and actuators to control or counter-act vibrations. One example of this type of prior art is U. S. 
Patent No. 4,490,841 which calculates Fourier transforms when in operation. 
3$ One important point to note about the above-listed approaches to active vibration control is that, although vibration 
and noise cancellation is a "control" problem, the above listed methods attempt to approach the subject of active vibration 
control from a "signal processing" and/or "filtering" view point. Thus, these techniques may not have a guaranteed 
stability. 

There is a need for an active vibration system that can be adapted quickly on-line to compensate for vibration, due 
40 to both random and repetitive excitations, in a system and which does not have the above-described drawbacks of the 
prior art. Further, there is a need for an active vibration control system which approaches the vibration control application 
from a "control" problem point of view. 

An improved active vibration system for controlling vibrations at a structure resulting from at least one excitation 
force acting upon the structure is provided. In at least one embodiment of the present invention the active vibration 
4S control system comprises at least one sensor which produces a sensor output and at least one actuator located at the 
structure for imparting a reaction force to the structure. Additionally, a controller Is connected between the at least one 
sensor and the at least one actuator, the controller Including a system identifier for receiving the output from the at least 
one sensor and deriving a relationship representative of a relationship between the output of the sensor and the appli- 
cation of the reaction force of the at least one actuator; and an optimal controller connected to the system identifier to 
BO receive the derived relationship and for developing control driving signals from the relationship for driving the at least 
one actuator. 

In one particular embodiment, the system identifier includes a neural network for learning the dynamtes of the struc- 
ture and for providing output signals that followthe state variables of the stmcture. In a first example, the optimal controller 
used involves Kalman filtering. 

ss In another particular embodiment, the optimal controller receives a state matrix A, an input matrix B and a state 

vector X from the system Identifier and uses Modified Generalized Predictive Control (MGPC) to derive said control 
driving signals for said at least one actuator. In a first particular example using MGPC in connection with the optimal 
controller, the system identifier uses least means square filtering to solve for the A and B matrices and the state vector 
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X. In another particular example, using MGPC in connection with tine optimal controller, the system identifier uses re- 
cursive least means square filtering to solve for the A and B matrices and the state vector X. In a further particular 
example, using MGPC in connection with the optimal controller, the system identifier uses a neural networl< to solve for 
the A and B matrices and the state vector X. 
5 Various embodiments of the present invention will now be described with reference to the accompanying drawings, 

in which:- 

FIGS. 1-4 show examples of systems wherein vibrations from random and/or repetitive excitations are acting upon 
a structure, those systems Incorporating the vibration control system of the present Invention; 

10 t 

FIG. 5 is a blocic diagram of a first prefen'ed embodiment of the present invention; 

FIG, 6 shows a block diagram of a vibration control system including a microprocessor, which system is suitable for 
use with the present invention; 

IS 

FIG. 7 is a block diagram of a first preferred embodiment of the present invention In which a neural network is used, 
on-line, to team the dynamics of the system; 

FIG. 8 Is a diagram of one possible implementation of a single neuron of the neural network shown in FIG. 7; 

so 

FIG. 9 is a diagram showing the interconnection of a plurality of the neurons of FIG. 8 to fomri the neural network 
of the present Invention; 

FIG, 10 is a schematic drawing of one possible hardware Implementation of the neural network of the present 

25 embodiment; 

FIG. 11 is a block diagram of the implementation of the algorithm of the neural network of the present invention 
used to solve for the matrices A and B and derive estimates of state variables; 

30 FIG. 1 2 Is a block diagram of the circuitry which may optionally be used to Implement the optimal controller of FIG. 
5 and/or the neural network and controller circuitry of FIG. 7 in one preferred Implementation of those embodiments, 

FIG. 13 is a block diagram which represents genericaily the ROM used in connection with the present invention, 
and the code stored therein which is used to implement the various embodiments of the present invention; 

3S 

FIGS. 1 3A - 1 3F are block diagrams representing genericaily the code stored In ROM of a microcontroller and which 
is used to implement various embodiments of the present invention; 

FIG. 1 4 Is a flow diagram showing the operation of the neural network used in connection with various embodiments 
40 of the present invention; 

FIG, 1 5 Is a flow diagram showing the operation of the modified generalized predictive controller for use with par- 
tk:uiar emlDodiments of the present inventbn. 

45 For the purposes of promoting an understanding of the principles of the invention, reference will now be made to 
the embodiment illustrated in the drawings and specific language will be used to describe the same, it will nevertheless 
be understood that no limitation of the scope of the invention is thereby Intended, such alterations and further modifica- 
tions in the Illustrated device, and such further applications of the principles of the invention as illustrated therein being 
contemplated as would nomnally occur to one skilled in the art to which the invention relates. 

50 Referring now to FIGS. 1-4, there are shown various examples of systems in which vibration or noise may occur. 
When a vibrating body, such as the vehicle 20 of FIG. 1 . the antenna 53 of FIG. 2. the engirie 63 of FIG. 3, or the motor 
80 of FIG. 4, is mounted on a support system by means of a mounting system, vibrations from the body are transmitted 
to the support system unless the mounting system is designed to either absorb the energy or counteract the vibrations. 
Additionally, certain systems may have further random vibrations associated therewith, such as vibration due to 

ss road inputs 30 of FIG. 1. The use of the particular systems shown in FIGS. 1-4 are meant to be for example only. 
Additional systems subject to vibrations include those Involving heavy machinery, sensitive computer and other electronic 
equipment, electric generators in an RV or on a ship or on board of an aircraft, fighter craft, or any civil structures and 
buildings exposed to earthquake and strong wind flows, stacks and mufflers, optical tables, medical equipment, en- 
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gine-body structure of a vehicle, etc., are all examples of systems with which the present Invention, as shown In the 
various embodiments herein, may be used. 

Referring more particularly to FIG. 1. there is shown a driver 10 seated inside a vehicle 20. In operation, various 
vibration promoting forces are at work on vehicle 20. For example, vehicle 20 Is subject to both road inputs 30, which 

s are random vibration sources, and engine torque pulses from engine 40, which cause repetitive vibrations. Vibrations 
resultant from the road inputs 30 and engine torque pulses from engine 40 can be relayed to the occupants of the 
passenger compartment. The resultant vibrations appear to the occupant(s) as either excess vibration and/or noise. 

Vehicle 20 additionally Includes an active suspension system Including, among otherthings such as shockabsbrbers 
48, the active engine mounts 42 which support the engine 40 on the vehicle body. 

TO Ideally, the optimal way to measure the forces acting on the vehicle 20 would be to measure road and engine 

excitations directly. Sensors so mounted are referred to as "direct feedfonward" sensors. However, It is impractical to 
put sensors on the tire or within the engine. As such, a vibration sensor 22, such as an accelerometer, a force transducer, 
a displacement sensor, etc., Is mounted on the steering column 26, and optionally, a microphone 28 is located near the 
driver's head area. Sensor 22 and microphone 28 are used to detect the efficacy or performance of the active vibration 

IS control system using feedback, and thus, are referred to herein as performance or feedback sensors. The performance 
sensors 22 and 28 relay a signal representative of the the amount of resultant vibrations detected in the passenger 
cabin 25 to the controller 43, which is, In turn, used to control the actuators of the active engine mounts 42. The actuators 
used with all embodiments of the present Invention may be of any type commonly known, some examples of which are 
electromagnetic, hydraulic, pneumatic, piezoelectric, electromechanical (forcers or shakers), electrohydraulic. shape 

20 memory alloy, etc. Performance sensor 22 produces a signal Indicative of the vibration sensed Inside the passenger 
cabin and relays that signal to the controller 43, Microphone 28 Is used to sense the amount of vibratfon detected In the 
passenger cabin 25 of the vehicle 20 which is manifested as noise. 

As noted herein, It would be useful to measure Incident excitations at th e excitation source, so as to provide advanced 
Information relating to the excitation (i.e. the magnitude and time of the resultant vibration). In many real systems, it is 

2S not possible to have a direct sensing of the vibration or noise source, or the source may not be identifiable as a localized 
excitation. However, since random excitations cannot be predicted from a simple feedback signal, then a direct source 
(feedfonrt^ard) measurement is needed. In order to resolve the problem of obtaining such a signal, this invention develops 
a control scheme that operates using a "pseudo'-feedfonward signal. Namely, using the example of FIG. 1 , if one cannot 
measure road excitations 30 on the tire 47. one can measure the excitation as far in advance of the structure (or per- 

30 fomiance point) as possible, i.e., use an accelerometer 49 on the axle near the tire and use that Information as a filtered 
(or "pseudo") version of the feedfonward signal to predict the required control action before vibration due to the road 
Input can reach the vehicle interior. Predictive measurement is possible, as the rate at which electronic transfer occurs 
Is much faster than the rate at which vibration or acoustic waves transfer through the applicable medium. 

As such. pseudo-feedfonA/ard sensors 44 and 49 have been provided to produce a filtered measurement of the 

3S' excitations close to their sources. More specifically, sensor 44 is mounted near engine 40, so as to produce and supply 
the controller with a signal Indicative of the repetitive excitations present near the engine 40. Likewise, sensor 49 is 
mounted on the axle 46, near the wheel 47 so as to provide a signal Indcative of the randomly occurring road inputs 30, 
acting upon the tire 47. As such, these pseudo-feedfonward sensors 44 and 49 have the ability to detect excitations prior 
to these excitations causing vibration and/or noise In the passenger cabin 25. The resultant signals measured by sensors 

40 44 and 49 are referred to as filtered or structurally filtered because they are not measured at the true excitation point, 
but rather at a point near the excitation source, and are thus filtered by the actual stnjcture connected between the 
sensor location and the excitation source. This is one example of a system using multiple pseudo-feedforward and 
feedback sensors. As will be seen in the examples herebelow, fewer or greater numbers of sensors and actuators may 
' be used in accordance with the present invention. 

45 Referring now to FIG. 2, there Is Illustrated an antenna 53, as a part of a stationary antenna system 50 which derives 
vibration control information from a performance (feedback) sensor 56. Between the antenna system base 51 and the 
antenna base 52, there are actuators 55 to impart vertical and horizontal forces in response to the control signals from 
controller 54. Controller 54 derives a response based upon input from performance sensor 56, which functions to detect 
the efficacy or performance of the active vibration control system using feedback. The antenna base 52 and Its associated 

60 antenna 53 can thus be stabilized effectively In a fashion similar to that set forth in connection with FIG. 1 . 

In FIG, 3, there Is illustrated an engine 63, as a part of a mobile engine system 60. Between the engine 63 and the 
motor mounts 62 on moveable frame 61 , there are actuators 65 to Impart vertical forces on opposite sides of the engine 
63 in response to the control signals from controller 64. Controller 64 derives control signals based upon information 
from the performance sensor 66 and the pseudo-feedforward sensor 67, located on the axle of the trailer 68 upon which 

ss the engine 63 is mounted. Performance sensor 64 f uncttons to detect the efficacy or performance ot the active vibration 
control system using feedback, whereas pseudo-feedfonA/ard sensor 67 sends a signal representative of a future vlbra- 
tbn due to road inputs or trailer vibration. The mobile engine system 60, thus has improved vibration control, in a similar 
fashion to that illustrated In connection with FIGS. 1 -2 and 4. 
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Referring now to FIG. 4. there is shown a fourth example of a system 70 In which vibration may occur. Consider 
sensitive electronic equipment 75, such as electronic navigation equipment, high precision machinery, or a computer, 
mounted to a support system 78, and having an active vibration control system in accordance with the present invention 
including actuators 76, on board a ship 71. In the system 70 of FIG. 4, vibrations from several different sources are 

6 acting on the ship 71 , and thus on the equipment 75 on board the ship 71 . A first source of vibration may be attributed 
to unpredictable (random) vibrations due to the ocean waves 72 striking the ship 71. As there can be no convenient 
direct measurement of these random vibrations available, a vibration sensor on the ship body, such as an accelerometer, 
can be used to pick up a filtered version of the wave effects. Additionally, the electric motor 80, which drives the rotor 
or propeller 79, provides a second, predictable source of vibration which acts upon the electronic equipment 75. 

10 In order to counter or control the amount of vibration Incident on the electronic equipment 75, the filtered or "pseudo" 

sensed signals representative of the effect of the waves 72 and the motor 80 on the ship 71 . are sent to controller 74 
from the sensors 73 and 81. Again, the signals produced by sensors 73 and 81 are ■pseudo'-feedfonward signals, as 
they are used to predict the future vibratory effect of the excltatk>n on the electrical equipment, before a resultant vibration 
is allowed to reach the electronic equipment 75. 

IS Additionally, a performance or feedback sensor 77, mounted at the support 78 and connected to the controller 74, 

relays to the controller 74 a feedback signal Indicative of the performance of the vibration control system. The signals 
relayed from sensors 73, 77 and 81 to the electronic controller 74 are used to predict the nature of the excitation and 
produce the required control signals that drive the active mounts 76 to protect the equipment 75 from the induced 
vibration. 

20 Altematively, sensor 81 may be omitted and repetitive vibrations such as are due to the motor 80 may be compen- 
sated for entirely using the feedback signal from sensor 77. Further, the signal from the sensor 77, may not only be used 
to predict and compensate for any repetitive excitations present in the system, it may also be used to eliminate any DC 
levels or leftover undesired vibrations and/or noise. 

Thus, the overall preferred system configuration of this Invention uses two signals: pseudo-feedfonward and feed- 

2S back, as shown in FIGS. 1, 3 and 4. However, it would be possible to implement the present invention using only a 
feedback sensor(s) as represented In FIG. 2. As stated above, the use of the examples shown In FIGS. 1 -4 is not meant 
to be limiting. A variety of systems subject to random and/or repetitive excitations would benefit by the use of any of the 
embodiments of the present invention. 

Referring now to FIG. 5. there is shown a block diagram of one implementation of the present invention, using a 

30 single pseudo-feedforward to detect random excitations and a single feedback sensor to report performance. This Is 
not meant to be limiting, as in various embodiments described herein whereas multiple sensors of each type may be 
used. I^ndom and repetitive vibrations, 502 and 504 respectively, add together, as represented by adder 505, and act 
on the structure 570. which, as noted herein, can be any structure subject to vibration and/or noise. In the example of 
FIG. 1 whereinasystemcomprisedofavehicle 20 was subjectedto random roadvibratbn and repetitive engine vibration, 

3S it was desirable to measure the random vibration as close to the excitation source as possible. This allowed for the 
prediction of a future vibration, prior to the actual vibration due to the excitation event being present, thus enabling the 
controller and actuators to counter the vibration before It reached the passenger cabin 25 of the vehicle 20. 

However, as noted in connection with that example, it Is Impractical to place the random excitation sensor directly 
at the source of random excitations, i.e. in the tire, so the sensor is located as close to the source as is practical. In the 

40 vehicle example, thefeedf onward sensor was placed on the axle, rather than the tire, as such the actual or direct excitation 
was not measured. Instead, a filtered excitation was sensed at the axle, which was the result of the actual excitation on 
the tire. As such, FIG. 5 shows the actual random excitation 502 being filtered by the structure of the system to which 
the pseudo-feedfonward sensor 520 is affixed, which Is represented by the structural filter 510. Pseudo-feedforward 
sensor 520 may be an accelerometer, a plezoceramic sensor, or other type of sensor as disclosed herein. Pseudo-feed- 

4S fonvard sensor 520 relays a signal representative of the filtered excitation to the system identifier 525. As used herein, 
the system Identifier may be either a system learning code or system learning device, as described herebelow, or any 
other system that Identifies the system characteristics. Additionally, a performance or feedback sensor 550, located on 
or near the stoicture 570, but away from the actuators, relays a signal which is representative of the sensed vibration 
present at the structure 570 to the system identifier 525. The object of the system identifier 525 is to derive a relationship 

so (i.e. a mapping function) between the signals obtained from the perfonmance and pseudo-feedfonward sensors 520 and 
550, and the application of the actuator inputs, control signals representative of which may be relayed to the system 
identifier 525 from the optimal controller 530 via feedback path 555. Altematlvely, infomiatlon relating to the actuator 
output may be obtained In other ways, for example, by direct sampling of the control driving signals using a sensor 
located at the actuator. The relationship between the application of the actuator outputs and the sensor inputs Is relayed 

ss to the optimal controller 530. 

Sensors 520 and 550, as discussed herein, may be multiple sensors or may each be a single sensor. Further, 
although not shown specifically in FIG. 5, a sensor may be placed near a source of repetitive excitations so as to provide 
a pseudo-feedfonward signal of that repetitive excitatton source (i.e., the sensor near the engine 40 In FIG. 1). 
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The optima! controllsr 530 uses a special optimal control algorithm and optimization to develop control driving signals 
for all of the actuators 560. Further, as described above, a feedback loop 555 provides a feedback signal from the optimal 
controller to the system identifier, for use by the system. The actuators 560 create the reaction or cancelling forces and 
apply them to the structure 570. Actuators 560 may be any known type of actuator, as described herein. FIG. 5 will be 
5 discussed in more detail herein In connection with various embodiments of the present Inventfon. 

Referring now to FIG. 6, there Is shown a block diagram of one possible Implementation of a vibration control system 

90 in accordance with the teachings of the present invention. FIG. 6 shows a system, such as the system described 
. above In connection with FIG. 5 or herebelow In connection with FIG.7, which may be used to Implement many of the 

embodiments of the present Invention described herein subject to any modifications or substitutions described herebebw 
10 In connection with each embodiment (i.e. the Implementation of the neural network using analog hardware rather than 
a microprocessor). The structure 99 represents the system structure which is subject to vibration from vibration source 
(s) 91 , and upon which the active vibration system of the present invention is designed to work. The vibration source(s) 

91 may be repetitive and/or random vibration sources, such as the road inputs 30 and engine 40 which act on the vehicle 
20 of FiG. 1. Sensor(s) 92a represent any pseudo-feedfonA/ard sensors that oiay be present in the system, (examples 

IS of which may be an accelerometer, a force transducer, a tachometer, a displacement sensor, etc.) and which are located 

near the excitation source. Performance 8ensor(s) 92b represent any sensors located on or near the structure 99 to 

provide feedback indicative of the efficacy of the vibration control system 90. 

In FIG. 6, both the system identifier (625 of FIG. 5) and the optimal controller (530 of FIG. 5) are implemented using 

a single controller 94, which receives the outputs of the sensors 92a and 92b. The outputs of sensors 92a and g2b are 
20 fed Into an analog signal converter/conditioning module such as A/D converter 93. If the signal from the sensor(8) 92a 

Is already provided in digital fomn, then the A/D converter 93 may be replaced by a filter for aliasing. 

The resulting digital signal is fed into the random access memory or RAM 95a associated with the controller 94. 

The stored data is processed by the processor or CPU and co-processor 96 of the system according to algorithms stored 

in the read only memory or ROM section 95 of the controller 94. These algorithms stored in the ROM section include 
2S the system teaming code, the optimization code, and the specialized optimal control code. These codes will be discussed 

with further particularity in connectfon with the descriptions of the various embodiments herebelow and in connection 

with FIGS. 13-13F. 

After proper data manipulation the GPU will generate appropriate control signals, which will be converted to analog 
signals using the D/A converter 97 and which will be sent to the actuator(s) 98. The actuator(s) 98 act upon the structure 

30 99 to counteract the vibrations caused by vibration source(s) 91 . 

Referring now to FIG. 7, there is shown a block diagram of a vibration control system 1 00 of a preferred embodiment 
of the present invention. In this first described preferred embodiment, the system identifier (525 of FIG. 5) comprises 
an artificial neural network, which is used to process sensory information from, feedback and pseudo-feedfonward sen- 
sors, so as to automatically learn on-line the dynamics of the structure, i.e., the vehicle 20, and to represent the learned 

35 dynamics in a state space form. 

One area of difference between the embodiments shown In FIG. 5 and that of FIG. 7, is that multiple pseudo-feed- 
forward sensors located near multiple excitation sources (which may be the sources of either random or repetitive vi- 
brations) may be used. This is because one advantage to the use of a neural network is that data from multiple input 
sources (i.e. pseudo-feedfonA^ard and/or multiple feedback sensors) may be processed using parallel distributed 

40 processing for fast learning and adaptation. Further, multiple outputs may be produced (MIMO) from the neural network. 
Although, this embodiment is shown using multiple sensors, in actuality, fewer sensors may be used. For example, one 
or more feedback sensors may be used in systems employing a neural network without pseudo-feedfonward sensors, 
or vice-versa. 

Vibration control system 1 00 includes pseudo-feedf onward sensors 1 44 and 1 49 (which may correspond to sensors 
45 44 and 49 of FIG. 1) located proximal to potential excitation sources. Additionally, perfonnance sensor(s) 124 (which 
may correspond to performance sensor 22 and microphone 28 of FIG. 1) are located on or near the structure 140 or 
the performance point of the system upon which the vibrations act. 

An artificial neural network 1 1 0 receives inputs from the pseudo-feedforward sensors 1 44 and 1 49 and the perform- 
ance sensor(s) 124. The neural network 110, acting to identify the system parameters, provides Inputs to an optimal 
60 controller 120 whbh analyzes the inputs and provides control signals to the actuators 114 to counteract the sensed 
excitatbns and vibrations acting upon the system 140. Feedback path 125 (similar to feedback path 555 in FIG. 5) 
provides a feedback loop from the output of the optinnal controller 120 to the input of the neural network 110, so as to 
provide additional Information to the neural network 110 to aid in learning the dynamics of the system. The adder block 
1 50 represents the summation of the vibrations due to the excitation sources and the control vibrations produced by the 
55 actuators to counteract vibration. 

The signals supplied to the neural network 110 from the sensors are used to "learn" the dynamics of the structure 
140 (which is. for example, vehicle, submarine, optteal table, industrial machinery, etc.) and to.provide output signals 
that exactly folbw state variables of the structure 140. These variables, including pertomiance measurements obtained 



6 



EPO 693 748 A2 

by the performance 8ensor(8) 124, are used by the neural network 110 to perform self adjustments. Knowledge of the 
various states enables the optimal controller 120 to provide a true optimum amount of required actuator force to coun- 
teract the incoming vibrations to the structure. 

A neural network, such as neural network 110, Is composed of a plurality of neurons. A biological neuron is believed 

s to be composed of a branching structure, comprising what are called dendrites, where the neuron Is believed to pick up 
signals from other neurons; a cell body, called the soma; a long transmission line-like structure called the axon, and 
brushlike structures at the tail end of the axon called synaptic buttons. One engineering interpretation of a biological 
neuron which may be used with the present embodiment of the invention Is shown in FIG. 6. 

Referring now specifically to FIG. 8, several input signals x^, X2, . . X|, . . Xn, are introduced to the input channels 

10 111 of the jth neuron 110a, Each individual neuron, of which neuron 110a is an example, takes each input with a diff eront 
weighting function 112, Wj,, and sends the weighted inputs through a summing amplifier 113, wherein the summation, 
U|, Is sent through a nonlinear tunctton [g] 114, typically a sigmoid function, to fomn the output signal, dj . The output 
channel 115 may be branched out and used to supply ah input signal to other neurons, as well as the output signal to 
the optimal controller 1 20 of FIG. 7, 

IS Mathematically, the processing of neuron 110a may be represented as 

wherein: 

20 S 

When a plurality of these neurons 1 1 0a are connected together, they form a neural network, such as neural network 1 1 0. 
The present preferred embodiment of the Invention uses an artificial neural network to represent vibrating structures 

2S and acoustic fields. Since each mode can be mathematically described by a second order differential equation, then 
two neurons are used for each mode in the structure or acoustic field. The first neuron is used to represent displacement 
and the second neuron is used to model the velocity of that mode. Thus for a structure having n modes of vibration. 2n 
neurons would be used. These neurons are connected in a manner known as the Hopfield neural-network, which has 
been described in, among others, papers entitled "Neural Computation of Decisions in Optimization Problems," Biological 

30 Cybernetics, Vol. 52, 1 985, pp. 1 41 -1 52, and "Simple 'Neural' Optimization Networks: An A/D Converter, Signal Decision 
Circuit, and a Linear Programming Circuit," IEEE Transactions on Circuits and Systems, Vol 33, No. 5, May 1986, pp. 
533-541 . Those articles are incorporated herein by reference. 

Referring now to FIG. 9 there Is shown one such structure for a neural network 1 1 0, based upon the Hopfield model, 
which may be used herewith in accordance with the present preferred embodiment of the invention. To model the fact 

3B that each neuron computes a nonlinear function of a host of inputs under the Influence of its own activation level, the 
electronic model provides a biasing current I] for each neuron. As such, FIG. 9 depicts n neurons U| - u^, each having 
an input bias signal l| - 1^, respectively. Additionally, each neuron Uj obtains an Input from the output from every other 
neuron U| in the network 110. as well as receiving a feedback input from itself. This feedback enabies the application of 
the Hopfield network to dynamic systems. Additionally, the output i)| of each neuron is provided to the optimal controller 

40 1 20. The equation of motion for each neuron's activity In continuous time for the network 1 1 0 of Fl G. 9, can be written as: 

^ ^>,=g(ii,) (1) 

^ * t WijVj + li (2) 

where Uj is the activation level at the i^h neuron, uj Is the output of the \^ neuron, and 1| is a bias term for the \^ neuron. 
W|j is the weighting for the connection between the ij^ and j^h neuron. 

Referring now to FIG. 10 there Is shown one possible way in which the the network of FIG. 9 may be constructed 
using electronic components. Each neuron may be made of a capacitor, resistor, and a nonlinear op-amp. The connection 
weights may be resistors of appropriate values. The implementation of the neural network 110 using the hardware 
described In connection with FIG. 10 Is further described in the above incorporated articles to Hopfield, et al., as well 
as others. Additionally, if the described neuron construction is used, an analog to digital (A/D) convector, not shown, 
^ should be used between the neural network 1 1 0 and the optimal controller 1 20 of FIG. 7. 

Although the neural network 110 may be implemented using analog hardware. It is preferrably Implemented by an 
algorithm stored in the ROM of a microprocessor, such as was discussed in connection with FIG. 6, and which will be 
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discussed more thoroughly In connection with FIG. 12. The use of a microprocessor with on-board ROM should not be 
limiting as either analog hardware or a microprocessor, microcontroller, digital signal processor, or a hybrid of analog 
and digital etc., may be used to implement all or part of the present embodiment. 

Further, in the present preferred embodiment, the microprocessor which implements the neural network 110 may 
additionally contain the optimal control code (as will be discussed In connection with FIGS. 1 3A, 1 3B and 1 3E). In which 
case, an A/D converter for receiving sensor inputs may be included in the circuitry leading into the microprocessor used 
to implement the neural network 110. 

In considering the structural dynamtos of the circuit of FIG. 1 0, consider the state space representation of the stmc- 
ture as: 

X = AX + By (2.3) 

Y = CX (2.4) 

where X Is a vector of the variables representing structural dynamics (e.g., modal characteristics) and V Is the excitation 
input. IVlatrices A and B contain structural parameters, e.g., natural frequencies and damping ratios. Y is a vector of 
structure output, representing miaasurements at performance points, C Is a matrix relating measurements and state 
variables. It should be noted that due to time varying characteristics of structures considered herein, the matrices A and 
B are assumed to be time varying. 

The neural network will be given vectors V and Y, where V is the outputs ot the sensors located near the excitation 
^ sources (i.e. 1 44 and 1 49 in FIG. 7) and Y Is the sensor outputs from the performance sensors (1 24 of FIG. 7). Additionally 
A(t) and B(t) must be identified in realtime. Thus, the neural network 110 (FIGS. 7-10) would automatically learn system 
dynamics using inputs from the sensors. As such, the present embodiment of the invention may be implemented on any 
vibration or acoustic system, because, given sensory Input, the neural network Is able to 'learn' the dynamics of the 
system. I.e. find the A and B matrices, and albw the system to automatically be tuned for the optimal active control of 
vibratiojtj or ng^ise cancellation. FIG. 11 shows a block diagram implementation of this invention's neural network 310, 
where A and B are estimates of the A and B matrices that the network is "learning". 

As stated above, the A matrix is a state matrix, while the B matrix is the input matrix. For every element in the A 
and B matrices, there will exist one neuron. For example, if the structure has three modes and two control actuators are 
used, then A will be represented by a 6 x 6 matrix and B will be represented by a 6 x 2 matrix. Thus, the artificial 
^ neural-network implemented will have a total of 48 neurons, 36 neurons for elements of the A matrix and 12 neurons 
for elements of the B matrix. The connection weights (W|j) and bias terms l| are derived on-line from the following equa- 
tfons. 

Wij(t) = -[Ti[*(£)<l. (£)e ^d^] (3) 



25 



40 



i.(t) - -[^f <&(£>[y(^)-ce x{0)]e d^l . (4) 

1 o 

n, = (A-KC)Qj+y,l (5) 

n, = (A-KG)n, + v,i (6) 

X = 0"^ e + exp (At)X(O) (7) 

_ A "a 

wherein A = A - KG, A = Estimate of A; and B s estimate of B. 

e%[ai,,a,2,-.a„„.b,,,bi2'- ^mJ (8) 



so 
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311 ai2 • • ' ^in 
A » 321 ^22 • • • 32n 



anl ^n2 • • • ^nn 

^11 ^12 • • • '^Im 
B « b2i ^22 • • • ^2m 



"^nl '^n2 • • • ^nm 



1 0 ... 0 

0 1 ... 0 s identity matrix 



again, wherein. W|j represents the connection weight (strength) between neurons \ and j, U| is the activation level of 
ss neuron j, -Oj is the output from neuron J, g( ) is a nonlinear, sigmoid, function of the neuron relating U and V. Further, m 
is the number of control actuators, n is the total number of state variables and l| is a bias term for a neuron. Additionally, 
in the above equations, A Is the state matrix. B is the Input matrix and C Is the output matrix. Additionally, r\ is the learning 
gain (0<ti>1), fi is the filter time constant, and, I is the identity matrix and Q, n are intermediate variables. Parameter 
\i Is chosen based on how fast the system dynamics are, because it represents how fast the previous data or measure- 
so ments should be forgotten and not to be Included in the calculation of the new estimates for A, B and X. FIG. 14 is a 
flow chart describing the operation of the neural network. 

Referring back to FIG. 7. In the case having the neural network implemented In software stored in the ROM of a 
microprocessor, as shown in FIG. 1 3A, the microprocessor would receive sensory data and the control output of the 
previous time-step and would perform the following: 

38 

(1 ) activate the neural network to use input/output data and learn the system dynamics. The results of this step are 
the calculation of the A and B matrices, and state vector X, as per FIGS. 11 and 14; 

(2) use the A, B, and C matrices to calculate the optimal control signals; and 

40 

(3) send the resulting control signals to the actuator(s) 11 4 by means of a digital to analog convertor 1 1 2. 

If the neural network is implemented in hardware, then the microprocessor only performs steps (2) and (3) above and 
an (A/D) converter Is provided between the neural network and the optimal controller to provide digitized inputs to the 
45 controller. 

More specifically, if the neural network 1 10 is implemented using a microprocessor, then, the microprocessor used 
to implement this neural network 110 and the optimal controller 120 operates as follows: 



(1 ) picks initial values for X(0), A(0) and B(0). as well as for a(0) and 11(0) (these may all In itialiy be chosen to be zero); 

(2) fomns the 6 vector; 

(3) selects the state variables and outputs according to the relationship Y = CX to identify the C matrix; 

(4) chooses the obsen/er gain K, which is typically five to 10 times the highest frequency considered for the structural 
analysis; 

(5) uses n and n to fomn the vector using equation (4.5); 
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(6) calculates the connection weights (Wij) and bias term l|, from equations (3) and (4); 

(7) reads signals from the sensors (1 44, 1 49 and 1 24 In FIG. 5) to obtain Uj and Yji 

(8) solves equations (1 ) and (2) based on the results obtained In steps 6 and 7 to find the neuron output; 

(9) sets the resulting n x n and m x n neuron outputs equal to the entries of the vector to find estimates for A and B; 

A A ... 

(10) calculates the Eigenvalues of A and computes K such that A-KC Eigenvalues are changed; 

(11) solves the first order differential equations (5) and (6) above to get Q and n; 

(1 2) finds the state variables, X, using equation (7) above. 

IS (13) repeats steps 5-1 2 above until the network reaches equilibrium; 

(14) stores A, BandX 

FIG 14 Is a flow diagram of the operation of the neural network 110 (FIG. 7) of the present embodiment in obtaining A. 

20 B and X. Since neural network has a unique, stable equilibrium state, the application of the method and apparatus of 
the present embodiment of the invention on a dynamic system, e.g., vibrating structure or acoustic field, will result in a 
true representation of the system. 

The objective of the neural network is to provide A and B matrices as well as the state vector X. Depending on the 
choice of states, the A matrix could directly provide natural frequencies and damping ratios of different modes. 

2S Referring back to FIG. 7, the information obtained by the neural network 110 as explained above, will be transferred 

to the optimal controller 120 logic portion of a microprocessor or microcontroller It should be noted that the neural 
network 110, can be used with a variety of control systems to develop control signals of different varieties for other 
control purposes as well as the vibration and noise control, as represented generally in FIG. 1 3A and more specifically 
in connection with FIGS. 1 3B and 1 3F described herebelow. The neural network of the present embodiment can addi- 

30 tionally be used for failure diagnosttes or estimating the state of any dynamic system of environment, as well as active 
vibration control. The optimal controller makes the foltowing assumptions; 

(1 ) that a maximum acceptable level of vibration and/or noise is provided (Yimax)) 

3S (2) that the maximum acceptable levels for the actuator outputs, e.g., forces or sound pressures (Vin,^^)* P^^' 
vided; 

Then, in one particular embodiment of the invention of FIG. 7, the optimal control code used in the optimal controller 
1 20, may be that known as Kalman filtering. Utilizing Kalman filtering, the optimal controller minimizes the following cost 
40 function: 

J » i l^tX^PX * it^Q it^at 

4S where P and Q are diagonal matrk:es that carry the terms 



and 



When J is minimized, then ratios of 



and 



Y^ 

'Imax 



^1 max 



Y, 



Y 

imax 
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V, 



V, 



imax 



5 



are minimized. 

To derive the optimal control input, V, the following equations are solved: 



r + ra + a'^r + rbq'"' b'^r-c'^pc = o 



10 



y = -Q'"" B*^ RX 

Therefore, the optimal controller 120 accomplishes the following steps: 



(10) 



(1 ) based on the and Ymax obtain Q and P matrices; 



(2) from the sensory Information regarding pertornrance points, obtain Y; 



IS 



(3) from the neural network receive A, B and X 



(4) solve the first order differential equation given in equation (9) above using any differential equatbn solver sub- 
routine, such as the known Runge-Kutta algorithm to find R; 



20 



(5) use equation (10) to derive optimal control signals to be sent to the actuators; and 

(6) repeat steps 2-6. 

^ As can be seen from the above, the embodiment of the present Invention shown in FIGS. 7-11 and 13A, 13B and 
13E includes the intelligence, using neural-networks, to learn the dynamics and system parameters of the vibrating 
structure or acoustic field automatically Additionally, as can be understood from the foregoing, the present embodiment 
provides for true optimal control of vibration and noise regardless of the number of sensors and actuators in the system. 
As such, the above described embodiment of the Invention can be applied to any vibration or acoustic system without 

^ a priori knowledge of the system. 

In accordance with another aspect of the present Invention there is provided an improved optimal controller which 
considers the limitatk>ns on the physical characteristics of actuator(s) in terms of the output level and the rate of change 
of the output. The presently described optical control system is described herelnbelow in connectbn with the model 
shown in FIG. 5. Note however, that it may be desirable to use the optinoal control system described herebelow in place 

^ of the Kalman filtering optimal controller 120 described in connection with the active vibration control system of FIG. 7. 
This Is because in typical applications of the optimal control using equations (9) and (10) as in the above described 
embodiment, the pseudo-f eedfonward signal is only used by the neural network to find A, B, and X, and not in the optimal 
controller 1 20. However, the present optimal control system further uses the pseudo-feedfonward signal In the architec- 
ture of the controller as well as in the computation of A, B, and X. 

^ The control technique used with a second embodiment of the present invention is an extension of the generalized 
predictive control presented by D. W. Clarke et al. in the paper •Generalized Predictive Control - Part 1 the Basic Algo- 
rithm," published in Automatica, Vol. 23, pp. 137-1 48,. 1987, that article being Incorporated herein by reference. The 
generalized predictive control (GPC) scheme has advantages over the use of other optimal control schemes. Whereas 
Kalman filtering is one optimal scheme for direct feedfonward and feedback signals, GPC or the GPC, as modified herein. 

^ (MGPC) is an optimal scheme for the general case which incorporates actuator output limitations and Its output rate 
limitations using feedback and pseudo-feedfonward signals. The Kalman filters and optimal regulators (LQR) (Equations 
(9) and (10) above) are most useful with a structure that has feedback only As such, the MGPC described herein is an 
improvement over the use of regular optimal control or Kalman filtering. 



As described above generally In connection with FIGS. 1 -4, both random and repetitive vibrations are known to act 



^ . on certain systems. Referring back to FIGS. 5, the present system has the advantage that It may be implemented using 
as few as two sensors, a single pseudo-feedforward sensor and a single performance sensor. Further, the pseudo-feed- 
fonward sensor(8) may be used to measure random and/or repetitive excitations. If more sensors are added, and the 
above described neural network embodiment were used in connection with an optimal controller using GPC or IS^GPC 
(FIG. 13F). the job of the controller would be simplified because A, B, and X estinnates may be provided faster. The 

^ actuators used may be of one or more of the types mentioned herein or known In the art. 

As described above, the system identifier 525 and the optimal controller 530 may both be Implemented using analog 
circuits, a microprocessor, microcontroller, a digital signal processor, or a hybrid of analog and digital processors in order 
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to provide flexibility in tlie use of mathematical expressions and in the number of input sensor signals and output actuator 

signals present. In connection with the present embodiment, a microprocessor 540 is used. FIG. 12 shows a more 

detailed view of one implementation of a controller, such as the controller 94 of FIG. 6, including a microprocessor 540, 

which may be used with the present Invention. The microprocessor 540.1s used to process sensor inputs according to 
s the code stored In the ROM 538 (FIG. 13). As noted above, a microprocessor 540 of the type shown in FIG. 12 could 

additionally be used to implement the neural network 1 1 0 and optimal controller 1 20 of the first embodiment (FIG. 1 3B), 

as well as with other optimal controllers (FIGS. 1 3A and 1 3B). 

Signals from the sensors 520 and 550 of FIG. 5 are first passed through a low pass filter 532 (FIG. 12) to eliminate 

high frequency noise. Then the signal is sent through an analog to digital (A/D) converter 534. Preferably, a 12-blt A/D 
10 converter is used in the present embodiment, however, other A/D converters may be used. It the neural network of FIG. 

7 Is implemented using a controller as In FIG. 12 then additional sensor Inputs may be provided. 

The signal from the A/D converter 534 is relayed to the microprocessor 540. As with many standard microprocessors. 

the microprocessor 540 includes a central processing unit or CPU 536, Including a clock, a math co-processor 537. 

on-board RAM 539 and ROM 538, which may take the form of EPROM or EEPROM. Further, communication (I/O) ports 
IS 541 and 542 are provided. As depteted in FIG. 13, codes stored in ROM 538 of the microprocessor 540 may provide 

for the system learning 543, optimization 544, and the optimal control algorithm 545. which in FIGS. 13C-13F is the 

modified generalized predictive control (MGPC) algorithm described herebelow. ROM 538 is accessed through interrupt 

signals from the CPU. 

The system leaming code 543 stored In the ROM 538, may be of many types. For example, in FIGS. 1 3A, 1 3B and 
20 13F, the system leaming code takes the form of the neural network described above. If the modified generalized predictive 
control system of the present embodiment is used as the optimal control in connection with a neural network, whether 
Implemented in code (FIG. 13F) or with hardware, the following equations provide the H(S) transfer function in the 
continuous domain and the IH(Z) transfer function in the discrete domain: 



2S 



30 



45 



BS 



H(S) = C(SI-A)'''b 
H(Z) = C(ZI-A)*"'b 



When a non-neural network system leaming code or device is used In connection with the MGPC. (i.e. FIGS. 1 3C 
- 1 3E) the following relationship is assumed between these points: 

b„+b-z'V...+b 

H(Z) = ^-1:^ ^— 

l+a^z +...+a„.,z 

where n is the order of the system and can be determined a priori by an experimental modal analysis of the particular 

35 system to which the Invention is applied. 

In order to guarantee consistency in the system learning, the controller sends a known signal, V(t), of high modal 
density (e.g. a dynamically rich signal) to the actuators, then simultaneously, sensor outputs (y(t)) are obtained and 
synchronized Input/output data streams would be stored in the RAM 539, after which, the system learning code would 
be called. If the system leaming code is not the neural network of the first embodiment (FIGS. 1 3c-1 3E) then a recursive 

40 fonnulation may be used. This recursive type system learning code, could use the Input/output data vectors to Identify 
afs and bj's of the system transfer function. The code contains the following two 
vectors and estimate of the output. y(t): 



A-f A A A A ^ .A, 

G =[3^,8^,32,---. a„..,.bT.b2.-- -.b J 
0^(t-1)=:[.y(t-1). -y(t-2),- - -y(t-n), V(t). V(t-1).- - -,V(t-n)] 



y(t) =G *<D 

where q is the estimate of unknown transfer function coefficients. 
BO Depending on the system complexity, two possible system learning algoritms useful with the MGPC control code 
may be either a least-mean-squares (LMS) or recursive-least-mean-squares (RLMS) algorithm (FIGS. 13D and 13E). 
Equations for these computations are set forth as folljjws;^ ^ 

LMS: ^(t)=ft(t-l)+27(y(t)-yp 
RLMS: e(t)=e(t-1)+K(t)[y(t)-y(t)] 



K{t) = P(klMdi 

Uo'^(t-1)P(t-1)«(t-1) 
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P(t)=:(1/X)[|.K(t)<D^(t-1)]P(M) 

where y is the identification gain (0<7<1 ). X is the data forgetting factor, and I Is an Identity matrix. After some transient 
period, q(t) and q(t-1) become very close meaning the system has been Identified, and the final values of a|'s and b|'s 
s are stored 

A least-mean-squares algorithm (FIG. 13D) is typically used for low order systems, whereas the recur- 
sive-least-mean-squares (FIG. 1 3E) algorithm is typically used for higher order stable systems. A neural network system 
as described herein may be used for ail types of systems adding intelligence to the system to guarantee numerical 
stability and unique solutions. 

10 The optimal control code and the optimization codes 544, are part of the optimal - predictive control algorithm of 

the present embodiment. The optimal predictive control algorithm of the present invention is a modification of the gen- 
eralized predictive control theory, in which that theory is made applicable to active vibratton and noise control problems 
using both pseudo-feedfonA/ard and feedback signals. 

Referring back to FIG. 5, the system 570. which Is exposed to random and repetitive excitations as represented by 

IS the adder 505, and additionally to control inputs from the actuators 560 can be mathematically represented as: 

A(Z-')y(t) = B(2"')V(t-1)+D(Z"^)7(t-1)+^ip5(t) 

where y(t) Is the performance signal from the performance sensor 550, V(t) is the active control input from the actuator 
20 (s) 560, 7(t) Is the random + repetitive signal once they have gone through the D(Z-'')/A(Z-^) filter to form the pseu- 
do-feedforward measurement, ^(t) Is the portion of the excitation that cannot be measured, and A is an integral control 
action forced on the control system to help eliminate the effects of the unmeasurabie random excitation, 4(t}. 
The MGPC optimal control signal V(t), is derived based on minimization of the following cost function: 



2S 



3S 



40 



so 



ss 



J-E {£ [y(t+j)-w(t+j>]2 + fn[AV(t+j-l)]2} (11) 

jaaO j=l 



Namely, It is the function of the control signal sent to the actuators to keep the performance point's vibration (or noise 
level) as close as possible to a defined threshold w(t), while minimizing the total Input effort, V(t). Furthermore, the 
30 parameter Fj is used to develop a desired trade-off between the control effort V(t) and the resulting performance Y(t). 

In selecting r, one is able to use W ^ , where Vj^gx is the maximum allowable control input. This aids in preventing 

max' 

the saturation of the actuator(s) 560. 

By considering the Diophantlne identity: 

C(Z'^) = EjAA+Z"'Fj 

The performance prediction y(t) is given by: 

y = G V + f 
where 

V - CAV(t), AV(t+l), . AV(t+n-l>]T 

^ where: 

G = E^B/C 



f=E|DA7(t+j-1)/C + F|y(t)/C 

y = [y(t). y(t+i).---,y(t+k)f 
k is the prediction horizon and the optimal control signal is 

V » (GTG+ri)-l GT(w-f) 
Note that v provides the true minimum for the cost function because: 
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a2j/av2 = (gTq + rii)T > 0 

which guarantees asymptotic stability of the controlled system. 

5 Theje are several design parameters that are selected to provide a fast and stable response, while satisfying the 

rate constraint of the control input, V(t). These parameters are the sampling rate, the prediction horizon for output (k) 
and for input (m), and the control penalty factor (r). In the preferred embodiment of the present invention, the eannpling 
rate has been chosen to be 1/1 0th of the settling time. This provides for a reasonable trade off between accurate rep- 
resentation of the Impulse response and the minimum computation time. 

10 The output prediction horizon (k) is based on the physical location of the pseudo-feedfonward sensor and the per- 

formance/feedback sensor. The selection of r and m are correlated such that as m is raised, the control signal becomes 
more activa As F is raised the control signal becomes less active. Thus m can be used to get close to the desired degree 
of control activrty, whereas r can be used for fine tuning. The controller 530 of the present embodiment of the invention 
uses the following rules to select m and y. 

IS 

(1) initially setm = r=0 

(2) raise m until AV(t) = rate constraint for 50% of the time 

20 (3) raise r until A V(t) = rate constraint for 33% of the time. 

For example, if the actuator is a shape memory alloy that can change its output (control force) at a rate of q Ibf/sec, 
then first, nri is increased until the actual and obsen/ed change of the output is 50% of g, that fixes the value of m. Then 
r Is raised until the observed changes are 33% of q. That determines the value of F. Thus there will be a short period 

25 of tuning of m and F during the initial stage of the design. Thus the above rules are stored in the ROM section of the 
microprocessor and as the system changes with time, new values of F and m are calculated automatically. 

In summary, the present invention relates to a vibration control system which may include a system identifier (system 
learning code or device) and an optimal controller. There is described herein one particular system Identifier which 
utilizes a neural network, which in one embodiment, is implemented using hardware. In another embodiment of the 

30 present invention, a neural network Is Implemented In software stored for use by a microprocessor or controller. The 
optimal controller utilized to process the outputs of the neural network was described in one particular embodiment as 
a Kalman filter algorithm stored in memory of a microprocessor or microcontroller. This Is not meant to be limiting as 
other optimal controllers may be used with the neural network described herein. 

Further, there is described herein an improved optimal controller which uses a OKxiification of generalized predictive 

35 control (MGPC) to process inputs provided to It by a system Identifier. In one particular embodiment the MGPC optimal 
controller was used to process inputs provided from a neural network, as described herein. In another embodiment, the 
MGPC optimal controller was used in connection with system learning code comprising a least-mean-squares compu- 
tation. In yet another embodiment, the MGPC optimal controller was described as being used in connection with system 
learning code comprising a recursive least-mean-squares computatbn. This is not meant to be limiting, as other types 

^ of system learning codes or devices may be used in connection with the MGPC scheme described herein. 

While the Invention has been Illustrated and described in detail in the drawings and foregoing description, the same 
is to be considered as illustrative and not restrictive in character. White there are shown in the drawings individual sensors 
and actuators, it should be understood that this is a schematic representation which can reflect a larger number of 
sensors and actuators. It is desirable to have at least three actuators and three sensors associated with the systems of 

4S the present invention. It is more preferred to have systems, for example, six sensors and three actuators. Having more 
sensors and actuators allows for compensation of additional modes of movement including rotational and translational 
motion about three axes. It is being understood that only the preferred embodiment has been shown and described and 
that all changes and modifications that come within the scope of the Invention are desired to be protected. 

so 

Claims 

1 . An active vibration control system for controlling vibrations at a structure resulting from at least one excitation force 
acting upon the structure, comprising: 
ss at least one actuator located at the structure for imparting a reaction force to the structure; 

at least one sensor located away from said actuator, said at least one sensor producing a sensor output; 
a controller connected between said at least one sensor and said at least one actuator, said controller Including: 
a system identifier for receiving said sensor output from said at least one sensor and deriving a relation- 
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ship between said sensor output and said reaction force imparted to the structure by said at least one actuator; and 
an optimal controller connected to said system identifier to receive said relationship and for developing 
control driving signals from said relationship for driving said at least one actuator; 

said system identifier including a neural network for learning the dynamics of the structure and for providing 
s ' output signals that follow the state variables of the structure. 

2. The active vibration control system of claim 1 . wherein said neural network comprises a plurality of artificial neurons, 
each neuron receiving a weighted input from every other neuron, as well as a weighted input biasing current related 
to the the output of said at least one sensor, said neural network providing an output from each neuron to said 

10 optimal controller. 

3. The active vibration control system of claim 2. which includes a plurality of sensors. 



4. The active vibration control system of claim 3, wherein said plurality of sensors includes at least one performance 
IS sensor located at said structure for providing an output indicative of the vibration sensed at said structure. 

5. The active vibration control system of claim 4, wherein said plurality of sensors additionally^ includes at least one 
pseudorfeedfonvard sensQ^bcated near the at least one excitation source. 

20 6. The active vibration control system of claim 5, wherein each said neuron additionally includes a feedback input from 
itself. 

7. The active vibration control system of claim 6, wherein each the weighting of said weighted inputs to said neuron 
is different for at least some of said inputs. 
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so 



8. The active vibration control system of claim 7, wherein said feedback input is also weighted. 

9. The active vibration control system of claim 7, wherein the output of each neuron ] is 

'Oj = g(Uj) 

wherein: 



n 

Uj = Z Wji^U). 

35 \ 

wherein g(} is a nonlinear function, \N^^ is the weighting function for the connection between the Jth and kth neuron, 
and is an input to the neuron. 

10. The active vibration control system of claim 9, wherein the output of said neural network are at least a state matrix 
40 A, an input matrix B, and a state vector X. 

11 . The active vibration control system of claim 1 0, wherein said optimal controller receives the state matrix A, the input 
matrix B and the state vector X from said neural network and uses Kalman filtering to derive saki control driving 
signals for said at least one actuator. 



1 2. The active vibration control system of claim 1 1 , wherein said optimal controller minimizes the following cost function: 



wherein P and Q are diagonal matrtees that carry the tenms 

1 



ss and 



I max 



*lmax 
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and, wherein eaid optimal controller additionally derives the optimal control input, SV, for said at least one actuator 
by solving: 

R + RA + a'^R + RBQ"'' b'^R - C*^ PC = 0 

S -1 T 

y = -Q^B^RX. 

1 3. The active vibration control system of claim 1 0, wherein said optimal controller receives the state matrix A, the input 
matrix B and the state vector X from said neural network and uses Modified Generalized Predictive Control to derive 

10 said control driving signals for said at least one actuator; 

14. The active vibration control system of claim 1 3, wherein said optimal controller minimizes the following cost function: 



IS 




where y(t) is the performance signal from the performance sensor, w{t) is a defined threshold, and V(t) is the input 
effort to said at least one actuator, Fj is a parameter used to represent a desired trade-off between the control effort 
V(t) and the resulting performance Y(t). 

so 

15. The active vibration control system of claim 7, wherein said neural network is ernbodied in software. 

16. The active vibration control apparatus of claim 7. including a plurality of actuators. 

25 17. A method for controlling vibrations at a structure resulting from at least one excitation force acting upon the structure, 
the structure being part of a vibration control system comprising at least one actuator located at the structure for 
imparting a reaction force to the structure, at least one sensor located away from the at least one actuator, the at 
least one sensor producing a sensor output, a controller connected between the at least one sensor and the at least 
one actuator, the controller including, a system identifier for receiving the output from the at least one sensor and 

30 deriving a relationship between the sensor output and the reaction force imparted to the structure by said at least 

one actuator and an optimal controller connected to the system identifier to receive the relationship and for devel- 
oping control driving signals from the relationship for driving the at least one actuator, wherein the system identifier 
includes a neural network for learning the dynamics of the stnjcture and for providing output signals that follow the 
state variables of the structure, comprising the steps of: 

35 

(a) using the neural network for learning system dynamics using input data from said at least one sensor and 
output characteristics and for deriving A, B and C matrices, and state vector X, wherein said A matrix is a slate 
matrix, said B matrix is an input matrix and said C matrix is an output matrix; 

40 (b) providing at least the A, B, and C matrices and state vector X to the optimal controller; 

(c) using at least the A and B matrices and state vector X to calculate the optimal control signals; and 

(c) sending the resulting optimal control signals to the at least one actuator. 

45 

18. An active vibration control system for controlling vibrations at a structure resulting from at least one excitation force 
acting upon the structure, comprising: 

at least one actuator located at the stmcture for imparting a reaction force to the structure; 
at least one sensor located away from said actuator, said at least one sensor producing a sensor output; 
so a controller connected between said at least one sensor and said at least one actuator, said controller including: 

a system identifier for receiving said sensor output from said at least one sensor and deriving a relation- 
ship between said sensor output and said reaction force imparted to the structure by said at least one actuator; and 
an optimal controller connected to said system identifier to receive said relationship and for developing 
control driving signals from said relationship for driving said at least one actuator; 
55 ■ said optimal controller using modified generalized predictive control incorporating actuator output limitations 

and output rate limitations in comblnatbn with the performance output and the pseudo-f eedfonivard output to derive 
said control driving signals for said at least one actuator. 
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19. The active vibration control system of claim 18. wherein said optimal controller receives a state matrix A, an input 
matrix B and a state vector X from said system identifier. 

20. The active vibration control system of claim 1 9, wherein said optimal controller minimizes the following cost function: 



where y(t) is the performance signal from the performance sensor, w(t) is a defined threshold, and V{t) is the input 
effort to said at least one actuator, Fj is a parameter used to represen- a desired trade-off between the control effort 
V(t) and the resulting performance Y(t). 

21. The active vibration control system of claim 20, wherein matrices A and B and state vector X are identified using a 
recursive method. 

22. The active vibration control system of claim 21, wherein said recursive method is a least means square (l^S) 
detenmination. 

23. The active vibration control system of claim 21 , wherein said recursive method is a recursive least means square 
(RLMS) determination. 

24. The active vibration control system of claim 21 wherein said system identifier includes a neural networlcfor leaming 
the dynamics of the structure and for providing output signals that follow the state variables of the structure. 



J-E {r [y(t+j)-w(t+j)]2 + Zri[AV(t+j-l)]2} 
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